----------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\Eric Zou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_figure2.log
  log type:  text
 opened on:  27 Jun 2022, 13:35:58

.         
.         use "$Rep_smokelabor/1_build/regdata/county_day.dta", clear

.         
.         tsset countyfip rfrnc_dt
       panel variable:  countyfip (strongly balanced)
        time variable:  rfrnc_dt, 01jan2007 to 31dec2019, but with gaps
                delta:  1 day

.         
.         * Wind dir. FEs
.         reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr g60_wdir) cluster(COUNTY10 rfrnc_dt) 
(dropped 38677 singleton observations)
(MWFE estimator converged in 13 iterations)

HDFE Linear regression                            Number of obs   =  3,968,857
Absorbing 3 HDFE groups                           F(  41,   1678) =      23.48
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.3770
                                                  Adj R-squared   =     0.2736
Number of clusters (COUNTY10) =      1,679        Within R-sq.    =     0.0416
Number of clusters (rfrnc_dt) =      4,269        Root MSE        =     4.9611

                  (Std. Err. adjusted for 1,679 clusters in COUNTY10 rfrnc_dt)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |
        F20. |   -.121052   .0792373    -1.53   0.127    -.2764663    .0343623
        F19. |   -.149058    .077071    -1.93   0.053    -.3002234    .0021073
        F18. |  -.2509384   .0809555    -3.10   0.002    -.4097229   -.0921539
        F17. |  -.2627857   .0831566    -3.16   0.002    -.4258873   -.0996841
        F16. |  -.1401913   .0825325    -1.70   0.090    -.3020687    .0216861
        F15. |   .0521488   .0853258     0.61   0.541    -.1152074    .2195049
        F14. |   -.034593   .0815367    -0.42   0.671    -.1945173    .1253314
        F13. |   -.033742   .0805736    -0.42   0.675    -.1917774    .1242935
        F12. |   .0245402   .0860235     0.29   0.775    -.1441844    .1932649
        F11. |   .0349867   .0820783     0.43   0.670        -.126    .1959733
        F10. |  -.0268984   .0842441    -0.32   0.750    -.1921329    .1383362
         F9. |   .0720917    .080375     0.90   0.370    -.0855541    .2297376
         F8. |   .0643738   .0809287     0.80   0.426     -.094358    .2231057
         F7. |  -.0559878   .0839678    -0.67   0.505    -.2206804    .1087049
         F6. |   .0687164   .0806959     0.85   0.395    -.0895587    .2269915
         F5. |   .1006141   .0846997     1.19   0.235     -.065514    .2667423
         F4. |   .1026689    .088677     1.16   0.247    -.0712603    .2765981
         F3. |   .2672466    .088611     3.02   0.003     .0934469    .4410463
         F2. |   .2683124   .0856667     3.13   0.002     .1002874    .4363373
         F1. |   .8987627   .0919531     9.77   0.000     .7184078    1.079118
         --. |   2.214289   .1036578    21.36   0.000     2.010977    2.417601
         L1. |   1.618362     .09738    16.62   0.000     1.427363    1.809361
         L2. |   .8481299   .0899668     9.43   0.000     .6716709    1.024589
         L3. |   .3399519   .0830314     4.09   0.000      .177096    .5028079
         L4. |    .339431   .0803612     4.22   0.000     .1818122    .4970498
         L5. |   .1947985    .084727     2.30   0.022     .0286168    .3609802
         L6. |  -.0080433   .0862084    -0.09   0.926    -.1771307    .1610441
         L7. |   .0087867     .08311     0.11   0.916    -.1542234    .1717968
         L8. |   .0794738   .0833521     0.95   0.340    -.0840113    .2429588
         L9. |   .1110324   .0848695     1.31   0.191    -.0554289    .2774937
        L10. |   .0363237   .0846749     0.43   0.668    -.1297558    .2024033
        L11. |    .036216   .0841099     0.43   0.667    -.1287552    .2011873
        L12. |    -.09392   .0834639    -1.13   0.261    -.2576243    .0697843
        L13. |   .0726878   .0814659     0.89   0.372    -.0870978    .2324733
        L14. |   .0812956   .0795565     1.02   0.307    -.0747448     .237336
        L15. |   .0811082   .0816452     0.99   0.321    -.0790289    .2412454
        L16. |  -.0171534    .085993    -0.20   0.842    -.1858183    .1515115
        L17. |  -.1529483   .0851182    -1.80   0.073    -.3198974    .0140008
        L18. |  -.0972914    .087989    -1.11   0.269    -.2698711    .0752883
        L19. |  -.1365937   .0808913    -1.69   0.091    -.2952522    .0220648
        L20. |  -.1312176   .0855015    -1.53   0.125    -.2989185    .0364833
             |
       _cons |   8.522859   .0385835   220.89   0.000     8.447182    8.598536
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
  county_doy |    564495      564495           0    *|
       st_yr |       637           0         637     |
    g60_wdir |         6           1           5     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         
.         * State x wind dir. FEs
.         reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr fe_wdir_st) cluster(COUNTY10 rfrnc_dt) 
(dropped 38677 singleton observations)
(MWFE estimator converged in 13 iterations)

HDFE Linear regression                            Number of obs   =  3,968,857
Absorbing 3 HDFE groups                           F(  41,   1678) =      23.37
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.3968
                                                  Adj R-squared   =     0.2966
Number of clusters (COUNTY10) =      1,679        Within R-sq.    =     0.0413
Number of clusters (rfrnc_dt) =      4,269        Root MSE        =     4.8819

                  (Std. Err. adjusted for 1,679 clusters in COUNTY10 rfrnc_dt)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |
        F20. |  -.1526101   .0764301    -2.00   0.046    -.3025184   -.0027018
        F19. |  -.1515763   .0740905    -2.05   0.041    -.2968959   -.0062568
        F18. |  -.2695156   .0782418    -3.44   0.001    -.4229775   -.1160538
        F17. |  -.2641965   .0802219    -3.29   0.001    -.4215421   -.1068509
        F16. |  -.1324182   .0793018    -1.67   0.095     -.287959    .0231226
        F15. |   .0553088   .0829301     0.67   0.505    -.1073486    .2179661
        F14. |  -.0379984   .0795308    -0.48   0.633    -.1939886    .1179917
        F13. |  -.0276159   .0785473    -0.35   0.725     -.181677    .1264451
        F12. |   .0374518   .0830577     0.45   0.652    -.1254558    .2003595
        F11. |   .0478926   .0800781     0.60   0.550    -.1091708    .2049561
        F10. |  -.0352481   .0813254    -0.43   0.665     -.194758    .1242619
         F9. |    .074693   .0775419     0.96   0.336     -.077396    .2267821
         F8. |   .0722918   .0786563     0.92   0.358    -.0819829    .2265665
         F7. |  -.0420898   .0822078    -0.51   0.609    -.2033304    .1191508
         F6. |   .0974807   .0779185     1.25   0.211    -.0553469    .2503084
         F5. |   .0993876   .0817831     1.22   0.224    -.0610201    .2597953
         F4. |   .0908916   .0862314     1.05   0.292    -.0782409    .2600242
         F3. |   .2488116   .0862793     2.88   0.004     .0795853     .418038
         F2. |   .2736606    .083637     3.27   0.001     .1096169    .4377044
         F1. |   .8851191   .0892926     9.91   0.000     .7099824    1.060256
         --. |    2.16414   .1009275    21.44   0.000     1.966183    2.362097
         L1. |    1.55238   .0956513    16.23   0.000     1.364771    1.739988
         L2. |   .8255589   .0876114     9.42   0.000     .6537198     .997398
         L3. |   .3591828   .0808116     4.44   0.000     .2006805     .517685
         L4. |   .3450962   .0781834     4.41   0.000     .1917489    .4984435
         L5. |   .1824861   .0819452     2.23   0.026     .0217606    .3432116
         L6. |  -.0028421   .0831215    -0.03   0.973    -.1658747    .1601906
         L7. |   .0219711    .080676     0.27   0.785    -.1362651    .1802072
         L8. |    .099467   .0807508     1.23   0.218    -.0589158    .2578499
         L9. |   .1090422   .0820892     1.33   0.184    -.0519659    .2700503
        L10. |   .0500205   .0820545     0.61   0.542    -.1109195    .2109605
        L11. |    .024934   .0825753     0.30   0.763    -.1370275    .1868954
        L12. |  -.0748677   .0807685    -0.93   0.354    -.2332854      .08355
        L13. |   .0625336   .0789622     0.79   0.429    -.0923411    .2174084
        L14. |   .0538653    .077181     0.70   0.485    -.0975158    .2052464
        L15. |   .0762283   .0796612     0.96   0.339    -.0800176    .2324741
        L16. |   .0012145   .0838576     0.01   0.988    -.1632619     .165691
        L17. |  -.1275475   .0830483    -1.54   0.125    -.2904366    .0353417
        L18. |  -.0947734   .0853359    -1.11   0.267    -.2621495    .0726027
        L19. |  -.1368957   .0780728    -1.75   0.080    -.2900259    .0162346
        L20. |  -.1309112   .0831116    -1.58   0.115    -.2939245     .032102
             |
       _cons |    8.52782   .0379188   224.90   0.000     8.453447    8.602193
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
  county_doy |    564495      564495           0    *|
       st_yr |       637           0         637     |
  fe_wdir_st |       294          49         245     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         
.         * County x wid dir. FEs
.         reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr fe_wdir_county) cluster(COUNTY10 rfrnc_dt) 
(dropped 38703 singleton observations)
(MWFE estimator converged in 21 iterations)

HDFE Linear regression                            Number of obs   =  3,968,831
Absorbing 3 HDFE groups                           F(  41,   1675) =      23.38
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.4064
                                                  Adj R-squared   =     0.3058
Number of clusters (COUNTY10) =      1,676        Within R-sq.    =     0.0416
Number of clusters (rfrnc_dt) =      4,269        Root MSE        =     4.8499

                  (Std. Err. adjusted for 1,676 clusters in COUNTY10 rfrnc_dt)
------------------------------------------------------------------------------
             |               Robust
        pm25 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    hms_deep |
        F20. |  -.1560764   .0761639    -2.05   0.041    -.3054628   -.0066899
        F19. |  -.1556051   .0737653    -2.11   0.035    -.3002871   -.0109232
        F18. |  -.2694386   .0779046    -3.46   0.001    -.4222393   -.1166379
        F17. |  -.2607498   .0798945    -3.26   0.001    -.4174533   -.1040462
        F16. |  -.1291906   .0789905    -1.64   0.102    -.2841211    .0257399
        F15. |    .059671   .0827308     0.72   0.471    -.1025957    .2219376
        F14. |  -.0360308   .0794675    -0.45   0.650     -.191897    .1198353
        F13. |  -.0274549   .0784345    -0.35   0.726    -.1812948     .126385
        F12. |   .0331346   .0829953     0.40   0.690    -.1296508    .1959201
        F11. |   .0483218   .0799734     0.60   0.546    -.1085365    .2051801
        F10. |  -.0354927   .0811586    -0.44   0.662    -.1946756    .1236902
         F9. |   .0817147   .0773835     1.06   0.291    -.0700639    .2334933
         F8. |   .0685196   .0785291     0.87   0.383    -.0855059    .2225451
         F7. |  -.0394766   .0820399    -0.48   0.630    -.2003881     .121435
         F6. |   .1006422   .0775967     1.30   0.195    -.0515545    .2528388
         F5. |     .09868   .0816055     1.21   0.227    -.0613795    .2587394
         F4. |   .0853677   .0859618     0.99   0.321    -.0832362    .2539716
         F3. |   .2497721   .0860735     2.90   0.004     .0809491    .4185951
         F2. |   .2694815   .0834358     3.23   0.001     .1058322    .4331309
         F1. |   .8797013   .0889333     9.89   0.000     .7052692    1.054133
         --. |   2.157388   .1005411    21.46   0.000     1.960189    2.354587
         L1. |   1.545086   .0952208    16.23   0.000     1.358322    1.731851
         L2. |   .8221966   .0873443     9.41   0.000     .6508812     .993512
         L3. |   .3612249   .0807233     4.47   0.000     .2028958    .5195541
         L4. |   .3519881   .0780977     4.51   0.000     .1988088    .5051674
         L5. |   .1825468    .081742     2.23   0.026     .0222195     .342874
         L6. |  -.0025758   .0829104    -0.03   0.975    -.1651947     .160043
         L7. |     .02143   .0805163     0.27   0.790     -.136493    .1793531
         L8. |   .0981795   .0804339     1.22   0.222     -.059582    .2559409
         L9. |    .112896   .0819816     1.38   0.169    -.0479011    .2736931
        L10. |   .0479047   .0818489     0.59   0.558    -.1126322    .2084416
        L11. |   .0231769    .082331     0.28   0.778    -.1383057    .1846594
        L12. |  -.0691892   .0804998    -0.86   0.390    -.2270801    .0887016
        L13. |   .0640836   .0788605     0.81   0.417     -.090592    .2187591
        L14. |    .050205   .0769339     0.65   0.514    -.1006916    .2011016
        L15. |   .0764091     .07956     0.96   0.337    -.0796384    .2324565
        L16. |     .00725   .0837707     0.09   0.931    -.1570562    .1715562
        L17. |  -.1284583   .0828946    -1.55   0.121    -.2910462    .0341296
        L18. |  -.0941764   .0850943    -1.11   0.269    -.2610788     .072726
        L19. |  -.1315598   .0778015    -1.69   0.091    -.2841582    .0210386
        L20. |  -.1329654    .082859    -1.60   0.109    -.2954836    .0295527
             |
       _cons |   8.527924   .0377655   225.81   0.000     8.453851    8.601996
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
     county_doy |    564482      564482           0    *|
          st_yr |       637           0         637     |
 fe_wdir_county |     10033       10033           0    *|
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         
.         ** plot
.         import excel using "$Rep_smokelabor/2_analysis/output_figures/appendix_figure2.xlsx", clear sheet("pm25_winds") first
(13 vars, 41 obs)

.         
.         tw      connected b all_wind  eday, col(blue "237 125 49") msize(0.8 1) msymbol(circle X) || ///
>                 line st_wind county_wind eday, col(black teal%80) lp(dash solid)  ///
>                 xtitle("Days since smoke exposure") ytitle("Change in PM2.5 (ug/m3)") ///
>                 xline(0, lcol(black) lp(dot) lw(0.6)) ///
>                 yline(0, lcol(black) lp(dot) lw(0.6)) ///
>                 ylab(, nogrid) ///
>                 legend(region(lcol(white)) order(1 "Baseline specification" 2 "Wind dir. FEs" 3 "State x wind dir. FEs" 4 "County x wind dir. FEs"
> )) ///
>                 graphregion(col(white))

.         gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure2.pdf", replace
(file C:\Users\Eric Zou\Dropbox\replicate_smokelabor/2_analysis/output_figures/appendix_figure2.pdf written in PDF format)

.         
. log close 
      name:  <unnamed>
       log:  C:\Users\Eric Zou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_figure2.log
  log type:  text
 closed on:  27 Jun 2022, 13:57:37
----------------------------------------------------------------------------------------------------------------------------------------------------
